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ABSTRACT 

Recent studies have indicated that the HCN-to-CO(J=l-0) and HCO+-to-HCN(J=l-0) ratios are 
significantly different between galaxies with AGN (active galactic nucleus) and SB (starburst) sig- 
natures. In order to study the molecular gas properties in active galaxies and search for differences 
between AGN and SB environments, we observed the HCN(J=l-0), (J=2-l), (J=3-2), HCO+(J=l- 
0) and HCO + (J=3-2) emission with the IRAM 30m in the centre of 12 nearby active galaxies which 
either exhibit nuclear SB and/or AGN signatures. Consistent with previous results, we find a sig- 
nificant difference of the HCN(J=2-l)-to-HCN(J=l-0), HCN(J=3-2)-to-HCN(J=l-0), HCO+(J=3 
2)-to-HCO + (J=l 0) and HCO + -to-HCN intensity ratios between the sources dominated by an AGN 
and those with an additional or pure central SB: the HCN, HCO + and HCO + -to-HCN intensity ra- 
tios tend to be higher in the galaxies of our sample with a central SB as opposed to the pure AGN 
cases which show rather low intensity ratios. Based on an LVG analysis of these data, i.e., assuming 
purely collisional excitation, the (average) molecular gas densities in the SB dominated sources of our 
sample seem to be systematically higher than in the AGN sources. The LVG analysis seems to further 
support systematically higher HCN and/or lower HCO + abundances as well as similar or higher gas 
temperatures in AGN compared to the SB sources of our sample. Also, we find that the HCN-to-CO 
ratios decrease with increasing rotational number J for the AGN while they stay mostly constant for 
the SB sources. 

Subject headings: galaxies: active — galaxies: ISM — radio lines: galaxies — galaxies: individual 
(NGC 1068, NGC 5194, NGC 4826, NGC 3627, NGC 4569, NGC 6951, NGC 6946, 
NGC 2140, M82, NGC6240, Mrk231, Arp220) 
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1. INTRODUCTION 

Activity in galaxies can be attributed to two main phe- 
nomena, highly active star formation, also called star- 
burst (SB), and mass accretion onto a supermassive 
black hole, often simply referred to as active galactic nu- 
cleus (AGN). Obviously, molecular gas plays not only 
a key role as fuel in the activity process but should 
also, in turn, be strongly affected by the activity. De- 
pending on the type, degree and evolutionary phase of 
the activity, different physical processes can be involved 
in changing the excitation conditions and chemical lay- 
out of the molecular gas, whether it is through strong 
ultra-violet (UV) or X-ray radiation fields or kinemat- 
ical processes such as galaxy interaction, large scale 
shocks, gas out- or inflow. Knowing the composition 
and characteristics of the molecular gas in active envi- 
ronments is thus essential for the understanding of the 
activity itself, its evolution and possible differences be- 
tween AGN and SB activity. Because of the differences 
in the radiation fields accompanying AGN and SB activ- 
ity, AGN are suspected to create excitation and chem- 
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ical conditions for the surrounding molecular gas sig- 
nificantly different from those in SB environments. In- 
deed, several recent studies, mainly based on molecular 
gas tracers such as CO, HCN, and HCO + (J=1-0), ap- 
pear to support this hypothesis: the HCN-to-CO(J=l- 
0) intensity ratios appear to be significantly higher 
and the HCO + -to-HCN(J=l-0) intensity ratios signifi- 
cantly lower in AGN (e.g., NGC 1068, NGC 6951, M51) 
than in SB (e.g . M82 NGC 6946) environment s (e.g 
Sternberg et all mffli iKohno et all Il999l 12001 

Kohno 2003, 2005). The difference in the intensity ratios 
can have various origins such as a) systematically differ- 
ent gas densities, b) systematically different gas temper- 
atures, c) different radiation fields (UV vs. X-rays) even- 
tually yieldin g different HCN, HCO+ and/ o r CO abun- 
danc e s (e.g.. [Tielens fc Hollenbach I 119851: iBlake et all 
19871 [Sternberg & Dakarno I I1995T fLepp fc Dalgarno I 
19961 : IMalonev et "aTlll996h . d) shocks, e) the evolution- 
ary stage of the activity, particularly important for star- 
burst, f) additional non-collisional excitation of the gas 
through IR pumping by UV/X-ray heated dust (e.g. 
iGarcia-Burillo et alj 120061: IWeifi et al J 120071 ). and g) su- 
pernova explosions (SNe), especially important for the 
HCO + excitation. 

In each case, the thermal and chemical structures of 
the gas should significantly differ between SB and AGN 
dominated regions. We thus carried out IRAM 30m ob- 
servations of three HCN and two HCO + transitions in 12 
nearby active galaxies (see Table [T]) to study the excita- 
tion conditions in SB and AGN dominated regions and 
their differences mainly as function of their gas densities, 
temperatures and molecular abundances. The sources in 
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this sample have been selected according to the following 
criteria: 1.) presence of either SB and/or AGN activity; 
2.) previously detected HCN(J= 1-0) emission; 3.) avail- 
able information on the CO emission; 4.) a declination 
above —20° so that they are observable from the IRAM 
30m telescope. 

2. OBSERVATIONS 

We observed HCN(J=l-0), HCN(J=2-1), HCN(J=3 
2), HCO+(J=1-0) and HCO+(J=3-2) in the centre 
(<30") of 12 nearby active galaxies 6 (TableQ] Fig. [Band 
[2]) with the IRAM 30m telescope at Pico Veleta (Spain) 
during January (HCN) and August 2006 (HCN+HCO+). 
In the January run, the AD set of SIS receivers was 
tuned in single side band mode to the redshifted fre- 
quencies of HCN(J=l-0) at 3 mm, HCN(J=2-1) at 2 mm 
and HCN(J=3-2) at 1 mm. In the August run, the AD 
set of SIS receivers were tuned in single side band mode 
to the redshifted frequencies of HCO + (J=1-0) at 3 mm, 
HCN(J=2-1) at 2 mm and HCN(J=3-2) & HCO+(J=3- 
2) at 1 mm. We used the 1 MHz backends with an effec- 
tive total bandwidth of 512 MHz at 3mm and the 4 MHz 
backends with an effective total bandwidth of 1024 MHz 
at 2 mm and 1mm. We spent ~l-4 hours on each target 
resulting in line detections with good to excellent signal- 
to-noise ratios (SNRs>5) for most of the sources. The 
atmospheric opacity at 225 GHz ranged between ~0.1- 
0.2 in ~80% of the time in January and between ~0.2- 
0.3 in August 2006. The (redshifted) HCN(J=2-1) line 
is still far enough in frequency from the atmospheric wa- 
ter absorption line at 183.3 GHz to be detectable with- 
out any contamination. Unfortunately, the same wa- 
ter absorption line prevents a reliable observation of the 
HCO + (J=2-l) transition, which is ~1 GHz closer to the 
water line. Also, the 2 mm receiver performance in the 
177.5-183 GHz window is severely reduced, yielding un- 
tolerably high system temperatures (>1500 K). We reg- 
ularly checked the pointing on a nearby planet and/or 
bright quasar resulting in a pointing accuracy within a 
few arcseconds (i.e., ~2-4"). 

3. DATA ANALYSIS 

Throughout the paper, the temperature scale is equiv- 
alent to T m b, i.e., main beam brightness temperature, 
which is defined as T m b=T*-F e ff/B ff . The beam (=B e g) 
and forward efficiencies (=F e gf) together with the beam 
sizes are given in Table [21 As beam filling effects are a 
crucial point for the analysis of our data, especially since 
we are mainly interested in the intensity ratios between 
different line transitions, all derived intensity ratios have 
been very carefully corrected using beam filling factors 
(see also Table [3]). We use the following definitions and 
relations: 
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Ju — 3,2; J i — J u 1; 
mol = HCN or HCO+ 



(1) 



with Rjj /iq being the line intensity ratio between 
the same molecule (i.e., HCN or HCO + ) at transition 
(J=J„— >Jr) and transition (J=l— >0), Ij ol / ; the intensity 

6 The HCO + emission from the three ULIRGs in our sample, 
namely Arp 220, Mrk 231 and NGC 6240, were taken from Gracia- 
Carpio et al. (2006). 



(see Equation[3]) of the molecule at transition (J=J U — >J;) 
and fj ul j, the filling factor (see Equation!!]). The HCO + - 
to-HCN intensity ratio at transition (J=J„— »J;) is de- 
fined as: 
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HCO+/HCN _ I? U C ?, + 
Ju.Ji ~~ If™ - ' 

Ju.J] 



J u — 3, 1; Ji — J u 1; 

(G=S + ) 



(2) 



I™ 1 ^ is the intensity of the molecular line in K kms 1 
at transition (J=J„— »J;) and defined as 

rmol _ P T ,mol(J=J u -J ' 



_ [■minunj=j u -j l j , 



Ju — 3, 2, 1; 

Ji = J u - 1; (3) 

mol = HCN or HCO+ 



fj u ,j, is the beam filling factor and defined as 

fju.Jl = 



(^"• J ') 2 
' Jl ) 2 +{< u 



Ju — 3, 2, 1; 

Jl = Ju - 1 



(4) 



with ^ u,Jl =size of the HCN(J=J U — Jj) emission re- 
gion (~size of the HCO + (J=J u — Jj) emission re- 
gion), #^ u,Ju ~ 1 = beamsize at the HCN(J=J U — Ji) 
(~HCO+(J=J u - Ji)) frequency (see TableQ; a is ei- 
ther 1 in case of a circular source or 0.5 in case of an 
elliptical source, i.e., one that fills the beam in one direc- 
tion (in the elliptical case we thus assume the minor axis 
of the emission as estimate for 8 S to derive the beam fill- 
ing factor while we use the Full Width at Half Maximum 
(FWHM) in the circular case). 

Some of the more distant sources in our sample are sig- 
nificantly smaller (<5") than the beamsizes resulting in 
smaller filling factors than for the more nearby galaxies 
(see Table 13]). The size of the HCN and HCO+ emission 
region for each galaxy has been estimated from inter- 
ferometric and/or single-dish HCN and/or HCO + maps 
where present or as a conservative upper limit from CO 
maps. Uncertainties of the order of 1-5" in the assumed 
source sizes translate into a <20% uncertainty in the in- 
tensity ratios. For most of the sources, we assume that 
the size of the HCN (HCO + ) emission region is simi- 
lar in all three (two) transitions, i.e., 9l° w dg 1 ~ df 2 . 
However, if this assumption were invalid, i.e., the size of 
the emission in the higher transitions is actually smaller 
by up to a factor of 50% than in the ground state, the 
estimated filling factor ratios between the different tran- 
sitions could be too high by up to a factor of 4. Thus, 
in the case of NGC 1068, we explicitly take into account 
that the HCN(J=3-2), HCN(J=2-1) and HCO+(J=3- 
2) emission are more compact (by a factor of ~1.5) than 
the HCN(J=l-0) and HCO+(J=1-0) emission as in- 
dicated by recent SMA observation of the HCN(J=3- 
2) and HCO+(J=3-2) emission in NGC 1068 (Krips 
et al. in prep.). As a conservative approach, we also 
use a lower size of the HCN(J=3-2), HCN(J=2-1) and 
HCO+(J=3-2) line emission for NGC 5194 whose in- 
tcrfcrometric maps also indicate a decreasing size of the 
emitti ng region with increasin g rotational number (Ta- 
ble M MiisuihitaiFin HSH . If we did not account 
for these differences in source sizes at different transi- 
tions, we might underestimate the line ratios for the AGN 
sources. Consequently, the so determined intensity ratios 
can be regarded as conservative upper limits for these two 
cases. 
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As we observed a large on the fly map (~2') for M82 in 
HCN and HCO + , we were able to average the emission 
in all line transitions over the same region resulting in 
identical filling factors. We thus assume a filling factor 
of 1 for M82. 

4. RESULTS 
4.1. The sample 

We detect all twelve galaxies in HCN(J=l-0) 
and HCO+(J=1-0), eleven in HCN(J=2-1), ten in 
HCN(J=3-2), and seven in HCO+(J=3-2) (Table G2 
Fig.[f]&[2}. The HCN intensity ratios are listed in Ta- 
ble Hand plotted in Fig. I3ll5l 

The diagrams (Fig. [3JSI) clearly indicate significant dif- 
ferences in the intensity ratios between the different ac- 
tivity types in our sample (see also discussion in the next 
sub-Sections): 

1. ) Rj°jJ is low (i.e., <0.4) in the 'pure' AGN sources 

of our sample and high in those with a dominant 
SB (i.e., >0.4), suggesting an increasing Rj^j* with 
an increasing SB contribution. The most extreme 
examples of R^ c £ are: NGC 1068 (low Rf™; 
AGN), NGC 6951 k Arp220 k NGC 6946 {moder- 
ate RH™; AGN+SB), and M82 (high Rf^; SB). 

2. ) Rfyio is low in thc 'P ure ' AGN sources (i.e., <0.3) 

but also in the composite (AGN+SB) sources, 
while it is high for the 'pure' SB sources (i.e., 
>0.4). In combination with R^/kp this creates 
so three different regions, separating the composite 
sources from the 'pure' SBs and AGN. The most ex- 
treme examples for each group are: NGC 1068 (low 

Rf^o & R32/10 ; AGN )> Ar P 220 & NGC 6951 

(moderate Rfffio k low RgS^,; AGN+SB), and 
M82 (high Rf 2 C /° + k R^S,; SB). 

3. ) R^ C< ^ ( / HCN is low to moderate in the 'pure' 

AGN and almost all composite sources (i.e., <1), 
while it is moderate to high for the 'pure' star- 
bursts (i.e., >1). The most extreme examples 

are: NGC 1068 (low R^° +/HCN ; AGN), Arp220 
(low R" C j ; +/HCN ; AGN+SB), NGC 6240 (high 
R 1°° +/HCN ; AGN + SB ), NGG 6946 & M82 ( hi 9 h 

■HCO+/HCN. 
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The apparent grouping of the two dominant activity 
types in our sample into different intensity ratios sup- 
ports fundamental differences between the excitation 
conditions of the two main activity types. It seems highly 
unlikely that biases in filling factors could lead to such a 
systematic trend. 

A comparison to CO data taken either from the litera- 
ture or from previous IRAM 30m observations reveals no 
similar separation effect in the CO line transition ratios 
for our sample (see Table HI . However, as a very inter- 
esting result, the HCN-to-CO luminosity ratios appear 
to decrease with increasing rotational number J in the 
AGN sources (Table [5J , while those of the SB sources 
remain more or less constant or even slightly increase. 



4.2. Individual sources 

In this Section, we will discuss some individual sources 
of our sample, each representing a good example of one 
of the ratio-extremes described in the previous Section. 

4.2.1. M82 - SB dominated galaxy 

M82 is the best testcase for a "pure" (evolved) SB in 
our sample. We mapped the entire central disk in M82 
with the IRAM 30m at 3 mm and obtained a number 
of discrete pointings along the disk at 2 mm and 1 mm. 
The HCN and HCO + data a re in good agreemen t with 
previous measurements (e.g.. iNguven et al] [T992V The 
gas disk in M82 is known to house two giant PDRs (e.g., 
iGarcia-Burillo et al.ll2002T ). prob ably also a central mas- 
sive black hole in for mation (e.g.. lMatsushita et alJ[2000t 
iPatruno et al. 120061) and a superbubble emerging fro m a 



past SNe fe.g. lWeifi et aUHOQalKronberg et al.l[l98lh . A 
different gas chemistry could hence be at play in the cen- 
tre than in the two PDRs. Thus, the HCN and HCO+ 
intensity ratios, discussed here, have been averaged over 
the entire map, and taken at the eastern PDR position 
and at the centre to allow for a better comparison. The 
30m observations indicate some variations especially of 

R?Vi0' R32/10 anc ^ ^32/10 between the position of the 

PDR and the nucleus (see Table HJ, while R^ ^ ^ HCN 
seems to be quite similar between the PDR and the cen- 
ter. The averaged ratios are very similar to the PDR 
ones indicating that PDR chemistry may dominate the 
overall emission in the disk. The PDR position in M82 

shows the highest Rf^ l/10 , Rf 2 C /° + and R* C ° +/HCN of 
all the sources in our sample (see filled red box in Fig. [3} 
[5]). The center in M82, however, seems to be, overall, 
similar to the other SB sources of our sample for most 
line ratios (compare next sub-Sections). 

4.2.2. NGC 6946 - SB dominated galaxy 

NGC 6946 is a local galaxy whose starburst activ- 
ity is assumed to be much younger than that in M82. 
NGC 6946 is not part of a gala xy-merger or -interaction , 
as it is the case for M82 (e.g.. iPisano k Wilcotsl 120001 ) . 
No signs of any significant PDR have yet been found in 
this galaxy and also large-scale, high-v elocity shocks do 
not s eem to play a major role yet (e.g.. ISchinnerer et al.l 
l2007t ). This difference to M82 in its SB properties may 
explain its location in the diagrams (Fig. [3][5]) with re- 
spect to M82, i.e., in the middle of the diagram, rep- 
resenting thus eventually a different type of starburst 
activity. As such, it might set tight constraints to our 
comparison between AGN and SB and underline the im- 
portance of the evolutionary stage of the starburst. 

4.2.3. NGC 1068 - AGN dominated galaxy 

NGC 1068 is the best example in our sample for hous- 
ing a pure AGN in a central radius of 1 kpc (=14") 7 . 
No strong evidence for any significant nuclear starburst 

7 Please note, that ^10-20% of the overall HCN emis sion in 
NGC 1068 is located in the spiral arms (Tacconi et al. 1994) which 
are possibly dominated by star formation. However, even for the 
largest beam-size of 29" the spiral arms should only marginally 
(~10%) contribute to the detected HCN emission in the 30m beam. 
The same is true for HCO + . 
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has been reported so far (e.g., MIR: Laurent et al. 
2000; NIR(PAH): Ima nishi 2002; Optical/Near- UV: Cid- 
Fernandes et al. 2001). iMarco k Brooks I (|2003ft estimate 
that a compact nuclear starburst would contribute less 
than 1% to the total IR luminosity. Thus, NGC 1068 
represents the best counter-part to M82 in terms of ac- 
tivity type. NGC 1068 appears to be always located in 
the opposite part of the diagrams in Fig. [3][5] with respect 
to M82, supporting the differences in the excitation con- 
ditions of the molecular gas sus pected between AG N and 
SB environments. Moreover, lUsero et al.l |2004) have 
discussed the possibility of NGC 1068 harbouring a giant 
XDR in its nucleus that is used to explain the suprisingly 
high HCN-to-CO and low HCO+-to-HCN luminosity ra- 
tios (see also Table [5]) . The potential prototypical XDR 
nature of the AGN in NGC 1068 classifies this source 
as an ideal counter-part to the PDR-dominated galaxy 
M82, in terms of effects of the radiation field onto the 
surrounding gas chemistry. 

4.2.4. NGC 5194 - AGN dominated galaxy 

Besides NGC 1068, the centre in NGC 5194 is most 
likely dominated by an AGN as well. This source is as- 
sumed to be in a post SB stage in which the massive star 
form a tion has already disappeared (e.g..lThronson et all 
1 199 It ISauvage et all 119961 : iGreenawalt et all I1998D : its 
nuclear activity is assumed to be caused by a low- 
luminosity active galactic nucleu s of LINER type (e.g., 
IHo. Filippenko. fcSar genii IT997T ). NGC 5194 is located 
close to NGC 1068 in all diagrams (Fig.JSJSl), substantiat- 
ing the differences between AGN and SBs in our sample. 
It is also one of the few sources for which a high HCN- 
to-CO but low HCO+-to-HCN(J=l-0)+ intensity ratio 
has been found (see Table [5]). 

4.2.5. NGC 6951 - AGN+SB galaxy 

The best example of a composite source in our sample 
likely is NGC 6951 (highlighted with a filled green tri- 
angle in Fig. [3][5|). It is known to house a prominent 
SB ring (e.g., optical: Marquez k Moles 1993; Woz- 
niak et al. 1995; Rozas, Beckman, k Knapen 1996; 
Gonzalez-Delgado et al. 1997; radio: Vila et al. 1990; 
Saikia et al. 1994) as well as a Seyfert type 2 /LINER 
nucleus in its central 20" (e.g., iBoer k Schulz 1 119931 : 
IHo. Filippenko. fc Sar genii 119971 ). iKohno et all (|1999l ) 
have mapped this source in HCN(J=l-0) with the 
Nobeyama array revealing that most of the HCN emis- 
sion is concentrated in the SB ring, similar to the CO 
emission. Probably due to missing sensitivity and angu- 
lar resolution (~ 3-4"), they fail to detect HCN(J=l-0) 
emission in the very centre (inner 2"). However, recent 
high angular resolution observations of HCN(J=l-0) in 
NGC 6951 carried out in the extended configuration of 
the IRAM PdBI show compa ct HCN emission in the cen- 
tre as well (jKrips et al.ll2007[ ) indicating a high HCN-to- 
CO(J=l 0) luminosity ratio of ~1. The interferomet- 
ric maps suggest that the HCN emission in the SB ring 
dominates the lines measured with the IRAM 30m tele- 
scope while the central HCN (and HCO + ) emission prob- 
ably contributes only marginally (less than 10%) to the 
observed brightness temperatures with the IRAM 30m. 
This explains the location of NGC 6951 in Fig. [3l close 
to the SB dominated sources in our sample. It has to 
be mentioned, though, that we might miss part (i.e., 



~30-40%)oftheHCN(J=3-2) and HCO+(J=3-2) emis- 
sion from the SB ring because of the lower beam size 
at these frequencies as recent SMA observation of the 
HCN(J=3-2) emission in NGC 6951 indicate (Krips et 
al., in prep.). However, this only slightly changes the 
position of NGC 6951 in Fig. MB The HCO+ emission 
separates it significantly from the evolved starburst M82, 
similar to NGC 6240 and Arp 220. This may be either 
a consequence of the evolutionary stage of the starburst 
or the additional existence of an AGN. 

4.2.6. Arp 220 - ULIRG (1 AGN + 1 SB nucleus) 

Arp 220 is the prototypical ULIRG and as such repre- 
sents well the higher-activity ULIRG population in our 
sample. A SB is ou tweighing the centre of Arp220 (e.g., 
iRisaliti et al.ll2006l and references therein) but an AGN 
may be present in at least on e of the two nuclei though 
this is still controversial (e.g.. ISanders Ifl9 88; Haa s et all 
120011 : llmanishi et al.ll2006l : iDownes k Eckartll2007t ). The 
central (i.e., 2-4") molecular gas emission in Arp 220 
is concentrated in two peaks, very similar thus in its 
gas morphology to NGC 1068. Despite the differences 
in origin for the gas morphology of these two galaxies , 
i.e., two nuclei in Arp 220 (e.g.. ISakamoto et al.1 119991) 
versus a (probably) w arped disk in NGC 1068 (e.g., 
iSchinnerer et al.l I2000T ) . the beam filling factor ratios in 
the 30m beam should be very similar for the two sources 
(when assuming similar source sizes in all transitions), 
making the differences between their spectra even more 
pronounced, especially in the HCN transitions (Fig. [TJ). 
Even if we account for smaller source sizes for the higher 
transitions in NGC 1068 (which actually leads to differ- 
ent filling factors between NGC 1068 and Arp 220; see 
Tabled, Rj u C £ /10 is still significantly lower in NGC 1068 

than in Arp 220. Interestingly, the situation for R^?o + 

and R^ ^ / HCN is much different: here the ratios of 
both sources are very similar to each other. This poten- 
tially differentiates Arp 220 significantly from M82 and 
the composite sources and, together with NGC 6240 (sec 
next sub-Section), may indicate that the ULIRGs have 
to be eventually handled as a different 'class' in our sam- 
ple, i.e., separate from the local starbursts and composite 
sources. Also, Arp 220 is located close to NGC 6946 in 
the Rj jj/io diagram (Fig. [3]) but farther away from it 

in those including HCO + . Its intensity ratios seem to 
be more like those of NGC 6951, which has a dominat- 
ing starb urst but also a cen tral weak AGN (see previous 
Section). lAalto et"ail (|2007j ) recently discussed the possi- 
bility of an XDR changing the molecular gas chemistry as 
an alternative to IR pumping to explain the high HNC- 
to-HCN(J=3-2) intensity ratios detected in Arp 220 by 
them. 

4.2.7. NGC 6240 - ULIRG (2 AGN) 

The two nuclei in NGC 6240 have been recently 
found to both harbor an AGN in addition to the 
very pronounced SB in this evolved merger (e.g., 
iKomossa et al.l I2003T) . Most of the molecular gas as 
well as d ust appears to be located between the two nu- 
clei (e.g.. [Bryant k Scoville 1 119991: iTacconi et all 119991: 
iNakanishi et al.1 l2005t llono et al.1 l2007h as opposed to 
Arp220 in which two gas disks may still be present (e.g., 
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ISakamoto et~aT1 Il999t ). Ilono et all (|2007f) also find evi- 
dence for a gas outflow/inflow that either is connected 
to starburst superwinds or outflows from the AGN. The 
existence of these and the two AGN in NGC 6240 may 
create exceptional and very extreme conditions for the 
molecular gas in the central region of NGC 6240. Sim- 
ilar to Arp 220, it is located close to NGC 6946 in the 
HCN diagram but closer to NGC 6951 in the diagram 

displaying R^o + (Fig. [J). However, in the R^ C ° +/HCN 

diagram (Fig. [5]), it seems to populate a very own region, 
quite separate from the rest, which may be linked to the 
extreme conditions in its centre. 

5. LVG SIMULATIONS 

We have run simulations of the excitation conditions 
for HCN and HCO + using the Large Velocity Gradi- 
ent (LVG) approximation in MIRIAD to connect the 
observed intensity ratios with physical parameters such 
as kinetic gas temperature, gas density and molecular 
abundances. Although we find differences in the val- 
ues obtained for HCO + between the LVG code used in 
MIRIAD and RADEX, we have concentrated our anal- 
ysis on the MIRIAD code. The differences between the 
two codes seem to be only present for gas regions with 
very high HCO + column densities (^10 15 cm -2 km -1 s) 
that lie, however outside the range studied in this pa- 
per. For lower column densities, RADEX and the LVG 
code in MIRIAD produce almost identical results. We 
think that the difference could be either a consequence 
of a different 'parameter sampling' between the codes or 
different assumptions that start to fail for higher column 
densities in one of the two codes. 

We assume a one component model for the LVG anal- 
ysis, i.e., all HCN as well as HCO + transitions originate 
from the same region underlying the same gas temper- 
ature and density. We carried out a reduced x 2 -test to 
constrain the above mentioned parameters; the x 2 ~test 

includes constraints based on Rj^/k)! Rfo/io ^ HCN and 

RHCO+ (TableHJ. 

We chose 4 exemplary sources in our sample, each rep- 
resenting one of the activity types and one of the ratio- 
groups: NGC 1068 as example of a pure AGN, NGC 6951 
as example of a composite source, M82 as example of a 
pure SB, and NGC 6240 as the example with the most 
extreme intensity ratios among the three ULIRGs in our 
sample. For the LVG analysis, we varied the gas tem- 
perature in a range of Tk=20-240 K using steps of 20 K, 
the H2 volume densities in a range of n(H2)= 10 1 ' 7 cm -3 , 
the HCN(J=l-0) column densities per velocity interval 
in a range of N(HCN(J=l-0))/dv~10 n " 19 cm- 2 km" 1 s 
and the HCO + -to-HCN abundance ratios in a range of 
[HCN]/[HCO+]=0.01-50. We also base the discussion 
and simulations on the following definitions and rela- 
tions: 

Z(HCN) = [HCN]/[H 2 ] = n(HCN)/n(H 2 ) (5) 

with Z(HCN) being the fractional abundance of HCN, 
[X] the molecular abundance of the molecule X and n(X) 
the volume density of X in units of cm -3 . We then have 



where N(HCN)/dv is the HCN column density per ve- 
locity interval in units of cm _2 /(kms _1 ), and dv/dr the 
"velocity gradient" in units of kms -1 /pc. 

Fig. [6] & [7] and Table [6] show the results of the LVG 
analysis with the lowest x 2 -values. Please note that 
we find two solutions respectively with similarly low 
X 2 -values for NGC 1068 and NGC 6951. The con- 
tours in the plots encircle regions with x 2 ^l an d are 
color-coded following the previous Figures (NGC 1068: 
blue; NGC 6951: green; M82: red; NGC 6240: yellow). 
The HCN abundance is additionally indicated in grey 
lines decreasing in thickness from 10 _5 (kms _1 /pc 2 ) -1 
to 10 _9 (kms _1 /pc 2 ) -1 in logarithmic steps of 1 based 
on equation 7. The standard (galactic) value for the HCN 
abundance found in galactic clouds is Z S (HCN)=2 x 10~ 8 . 
Giving a velocity range of ~100-500 kms -1 (compare 
Table [3]) for this sample and assuming typical sizes 
of —100-500 pc, we can constrain dv/dr to be in the 
range of ~0.2-5kms _1 /pc. This results in a range of Z a 
(HCN)/(dv/dr)~10~ 9 -10- 7 (kms-Vpc 2 ) -1 if assuming 
standard HCN abundances. However, for the AGN 
sources in our sample dv/dr is probably rather —1-5 
(large velocity width (—300 kms -1 ) and quite compact 
gas regions (<100 pc)) as opposed to —0.2-1 for the (lo- 
cal) SB(+AGN) dominated sources (covering a more ex- 
tended region of a few 100 pc than the AGN; e.g., M82 
or NGC 6951), i.e., 



Z s (HCN)/(dv/dr) AG N < 2 x 10~ 8 (km s" 1 pc 2 )" 1 



while 
Z s (HCN)/(dv/dr) SB 



> 2 x 10~ 8 (km s- 1 pc 2 )- 1 



(7) 



(8) 



N(HCN)/dv = Z(HCN)/(dv/dr) x n(H 2 ); 
with N(HCN) = n(HCN) x dr 



(6) 



if taking a standard Z S (HCN). 

The LVG analysis most strikingly constrains the molec- 
ular gas densities in the AGN and SB(+AGN) sources. 
While HCN and HCO+ emission in the SB domi- 
nated sources, such as M82 and NGC 6951, seem to 
emerge from regions with high H 2 densities in the 
range of n(H 2 )=(10 4 -10 6 - 5 )cm~ 3 , the HCN and HCO+ 
emission in sources of our sample with a pure AGN 
seem to be restricted to regions with gas densities of 
n(H 2 )<10 45 cm~ 3 . The LVG analysis also restricts the 
kinetic gas temperatures to be below <120 K for the 
SB sources while the pure AGN cases appear to have 
no upper limit for T^. The possibility of significantly 
higher temperatu res in AGN has b e en alr eady discussed 
for NGC 5194 bv lMatsushita et all (|1998l ) to explain the 
13 CO emission in this source. A high kinetic tempera- 
ture of T k >70 K has been also r eported for NGC 1068 
before by Sternbe rg et al.l (|1994h . We furthemore find 
tendentially larger HCN abundances in the AGN sources 
than in the SB sources which may even lie significantly 
above Z S (HCN) by a factor of -2-200. M82 thereby 
seems to denote a different extreme, i.e., having an HCN 
abundance that appears to be lower by at least an or- 
der of magnitude than Z S (HCN). We also find an ex- 
tremely low [HCN]/[HCO + ] abundance ratio of <1 that 
implies a very low HCN abundance and/or an addition- 
ally increased HCO+ abundance. The [HCN]/[HCO+] 
abundance ratios seem to be around 10 (or more conser- 
vatively: 1<[HCN]/[HCO+]<50) for the rest. 

The LVG analysis does not yield a very good solution 
for the HCO + emission in NGC 6240, which seems to 
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stand out a little bit with its intensity ratios, especially 
in R^/io ^ HCN , compared to the other two ULIRGs. The 
intensity ratios in Arp 220 and Mrk 231 are more like 
NGC 6951 in their ratios. This might indicate that a 
one-component model may be too simplistic in this case, 
whil e it reproduces wel l the values of the rest of the sam- 
ple. iGreve et~a l. (2006) present a two-phase LVG model 
for NGC 6240 (and Arp 220) that seems to fit nicely their 
data including also several transition of HCN and HCO + . 
However, our results on the gas densities and tempera- 
tures are not inconsistent with their findings, although 
our n(H2) seems a little bit lower than theirs; this may 
be though a consequence of our one-component model 
as their two phase model assumes a very dense plus a 
moderately dense gas phase which should result in an 
averaged and thus lower density in a one-phase model. 

The systematically increased HCN abundance 
in the AGN sources of our sample implies that a 
mass determination through HCN in AGN may 
significantly overesti mate the d e nse m olecular gas 
mass (compare also iKrips et al.l |2007( ) when using 
the same conversion factor as for SBs/ULIRGs. 
Based on the LVG analysis and adopting HCN 
abundances in the range of Z(HCN)/(dv/dr)w(0.1- 
10)xl0~ 7 pc/(kms- 1 ) with T k =40-240 K for AGN, 
we find brightness temperatures (T b ) in the range 
of -10-100 K and H 2 gas densities of n(H 2 )=10 20 - 
10 45 cm -3 . Taking the molecular mass (Mh 2 [M q ]) 
to HCN lum inosity ratio (Lhcn[(K km s" 1 pc)1) 
of (e.g., [Solomon. D ownes & Ra dford I Il992t 
ISolomon. Radford fc Downes 1119901 ): 

X(HCN) = M H2 /Lhcn = 2.1 x n(H 2 ) a5 /T b , (9) 

we obtain X(HCN)=10t7° M (K km s" 1 pc 2 )" 1 for 
the AGN sources. This conversion factor is ~2 
times smaller than that d e rived for ULIRGs by 
ISolomon. Downes fc Radford I 0-992) and IGreve et al.1 
(l200l of X(HCN)=20± 3 ° Mq (K km s" 1 pc 2 ) " 1 but 
agrees with that favored by iGao fe Solomon 1 |2004a) 
for their sample of nearby active galaxies and ULIRGs. 
An estimate of X(HCN) for the SB dominated sources 
in our sample yields a similar value to that of 
ISolomon. Downes fe Radford I (|1992D and IGreve et al.1 
(200^ _ 

By comparing the brightness temperatures with the 
obtained kinetic temperatures (for the J=l-0 lines), we 
can give a very crude estimate of the optical depth 
r, assuming r=-ln(l-Tb/T k ). This yields r«0.4-1.0 
for the SB dominated sources, r«0.3-1.4 for the AGN 
sources, r«0.3-1.8 for the composite sources and t«0.7- 
1.0 for NGC6240, indicating similar opacity ranges. 
However, the ratio between T b and T^ could be eas- 
ily biased if the emission is very clumpy, artificially 
lowering the T b /Tk ratio. As we probably average 
over many giant molecular cloud complexes, a clumpy 
structure may be ind eed s uspected in both SB galax- 
ies Ce.g.. IZhang et al.ll2001t lAlonso-Herrero et alj l2002t 
IWilson et al.ll2003tlGalliano et al.ll2005D and also for the 
gas/dust emission around AGN. 

6. DISCUSSION 
6.1. Chemical layout and excitation conditions 



We briefly described various physical processes in the 
Introduction that can lead to the observed differences 
in the line and transition ratios between AGN and SB 
dominated galaxies. As first scenario, we mentioned 
gas densities and temperature effects (cases (a) and (b) 
in the Introduction): higher gas densities and/or tem- 
peratures tend to increase populations of higher-J CO 
levels and may lead, for a given column density, to a 
reduced CO(J=1-0) line intensity (assuming that the 
gas densities are not that high to have a similar ef- 
fect on the HCN emission). Given the high HCN-to- 
CO(J=1-0) intensity ratios in AGN, this would conse- 
quently mean that the gas densities and/or temperatures 
in AGN must be higher than around SB activity. We can 
clearly exclude higher (average) gas densities in AGN, 
as our LVG analysis yields quite 'low' gas densities of 
n(H 2 )<10 4 5 cm -3 in the AGN sources. However, within 
the 30m beam, we potentially average over many giant 
molecular cloud complexes (GMCs) which could each ac- 
tually be denser than their average, if, for instance, the 
CMC number density is not very high and lower than 
that in the SB sources. Thus, one explanation for the 
'lower' average gas densities may be that the molecular 
gas traced by the HCN and HCO + emission does not 
split up into as many high density clumps/ GMCs in our 
AGN sources as it seems to be the case in SB environ- 
ments (compare IZhang et al.ll2001l:IAlonso-Herrero et al.l 
l2002HWilson et al.ll2003l: iGalliano et al.ll2005h . 

An increased HCN abundance in AGN compared to SB 
environments could equally explain the observed differ- 
ences in the line ratios. A chemical enhancement of HCN 
can be created in two different ways: either through far- 
ultraviolet radiation from O and B stars in young high- 
mass star forming reg i ons (e.g.. iTielens fc H ollcnbach 
[l98llBlake eTaII[l987t ISternberg fc Dakarno Hl995ir or 
through strong X-ray (i r )radiation from an A GN (e.g., 
iLepp fc Dalgarno 1 119961 : iMalonev et "all 119961 ). While 
UV radiation affects primarily the surfaces of gas clouds 
in the circumnuclcar regions (<lkpc), creating Photon 
Dissociation Regions (PDRs), X-rays penetrate deeply 
into the circumnuclear disk (CND), forming huge X- 
ray Dissociation Regions (XDRs). As a consequence 
of this volume versus surface effect, the X-ray radia- 
tion from the AGN might thereby produce higher HCN 
abundances relative to CO than the UV radiation of 
SBs which may explain the significantly higher HCN- 
to-CO(J=l-0) lum i nosity ratios found in AGN (e.g., 
Kotmo et all 120011: lUsero et alj 120041 : llmanishi etall 
20061 : Gracia-Carpio et al.l |2006[ ). This scenario seems 



to be consistent with the estimated abundance ratios in 
our sample. They appear to be higher in the AGN than 
in the SB dominated sources in our sample. 

Alternatively, the CO abundance might be smaller 
in AGN than around SB a ctivity either because 
of oxygen depletion (e.g. , ISternberg et al.1 119961 : 
iShalabiea fc Greenberg 1996) o r the influence o f X-rays 
from the AGN fe.g iMeiierink fc Spaans I [2005h . This 
could also result in the apparently higher HCN-to-CO 
ratios in AGN. The oxygen depletion has, however, been 
already ruled out for the XDR in NGC 1068 because it 
would produce different HCO+-to-HCN(J=l-0) ratios 
than observed (e.g., Usero et al. 2004; Kohno et al. 
2005). On the other hand, at the presence of a strong 
X-ray radiation field, CO dissociation may occur more 
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frequently predicting similar HCO + -to-HCN intensity 
ratios to those observed in NGC 1068 and NGC 5194 
and eventually also a de creased CO abundance (e.g. 
iMeiierink fc Spaans I I2005T) . This scenario can neither 
be supported nor discarded with our current data and 
has to be investigated further. 
Supernovae explosions (SNe), or more generally, ionisa- 



tion effects from cosmic rays (e.g., iDickinson et alj|1980t 



iWooten lUMH : iNguven et aDljofa MM et al.lll992f L are 

suspected to significantly increase the HCO + abundance 
while potentially decreasing the HCN abundance, yield- 
ing thus higher HCO + -to-HCN intensity ratios in evolved 
SBs (higher frequency of SNe) than in AGN. The most 
prominent examples for the role of the evolutionary 
stage of a SB are M82 and NGC 253; the starburst in 
NGC 253 is supposed to be in an evolutionary stage 
prior to that one in M82 and its molecula r gas appears 
to be mainly dominated by shocks (e.g., iMartin et al.l 
20061) rather than P DRs/SNe in contrast to M82 (e.g., 



Garcia-Burillo et al J 12002) . The HCO+-to-HCN inten- 



sity ratio is observed to be significantly lower in NGC 253 
than in the PDR of M82 (see Fig. 0§. In this context, 
the results obtained on the HCN and HCO + abundances 
in the evolved SB, PDR-dominated M82 are in excellent 
agreement with what is expected from theoretical predic- 
tions; M82 shows a very low [HCN]/[HCO + ] abundance 
ratio of <1. This may also explain the difference to the 
other SB sources in our sample who may all be in an ear- 
lier SB phase than M82, and emphasize the dependence 
of the excitation and chemistry on the evolutionary stage 
of the SB, whether reflected in more pronounced shocks 
or PDRs/SNe dominance. This is of particular impor- 
tance as it sets tight constraints on the comparison to 
AGN dominated sources. 

As a further alternative, dust heating, coupled to gas 
density and caused by the UV/X-ray radiation of the 
AGN, together with infrared (IR) radiative pumping 
could also be the origin of the stronger HCN emission 
in XDRs although, if true, a tighter correlation of the 
HCN luminosity with the IR luminosity than with the 
FIR luminosity would be expected as well as a strong 
correlation between the IR and X-ray lumin osity. Such 
correlations have not yet be en found (e.g., iLutz et all 
l200ilOao fc Solomon ll2004alrbT ) but X-ray absorption in 
compton thick sources and variability effects might lead 
to a large scatter in the data used so far, thus washing 
out possible correlations. However, IR pumping should 
similarly affect the HCN and HCO + emission, which 
both have vibrational modes at si milar wavelengths o f 
12-14 /im, as recently discussed by (jGuelin et a.1.1 l2007h . 
Thus, HCO + -to-HCN ratios close to unity would be ex- 
pected, in case IR pumping is significant, implying that 
non-collisional excitation of HCN and HCO + in sources 
with low HCO + -to-HCN ratios cannot be a dominant 
process. Also, the HCO + -to-HCN intensity ratios in the 
AGN sources even decrease with increasing J; this would 
not be expected if IR pumping were significant. How- 
ever, recent observations of the HNC emiss ion in some 
ULIRGs and a high re dshift quasar (e.g., lAalto etlrtl 
120071 IGuelin et ai1l2007h indicate that IR pumping may 
well be a significant factor in the excitation of at least 
HNC in sources with compact IR nuclei, but HNC (its 
vibrational state is at 24 /im) seems to be more easily ex- 
citable through IR pumping than HCN and HCO + (e.g., 



lAalto et al.ll2007t IGuelin et al.ll2007l ). HNC observations 
of our sample will be discussed in a future paper. 

We note that, although the statistics in our sample 
may be still small, especially for the AGN sources, the 
trend of different line ratios between AGN and SB dom- 
inated sources, seen in our data, is very pronounced and 
also further supported by literature data. Not only fall 
sources such as NGC 253 and NGC 4569 in the same 
line ratio regions as our SB sources, but the results on 
the AGN sources are also e ncouraged by rec ent results 
on the center of NGC 6951 (|Krips et al.ll2007l) as well as 
that on NGC 1097 and NGC 5033 (e.g. iKohno 1 120051 
iKohno. Nakanishi. fc Imanishil [2007). However, more 
AGN dominated galaxies in our sample would certainly 
be beneficial to substantiate our results. 

6.2. Comparison to theoretical results 

The lower gas densities in our AGN sources may actu- 
ally res olve a 'conflict' introduced by recent theoretical 
studi es (jMeijerink fc Spaans ll2005tlMeijerink et al.ll2006l . 
I2007T ) on XDR and PDR environments. While these au- 
thors claim that the HCO + -to-HCN intensity ratios are 
actually >1 in XDRs, we rather find low HCO+-to-HCN 
intensity ratios as opposed to their results. However, 
as emphasized in Mcijcrink et al. (2007), high HCO + - 
to-HCN intensity ratios are expected in high density re- 
gions, which does not seem to be the case for the 100 pc 
scale disks/regions in NGC 1068 and NGC 5194. For 
lower density regions, Meijerink et al (2007) indeed find 
low HCO + -to-HCN intensity ratios in agreement with 
our results. Also, their models predict higher surface 
temperatures in XDRs for strong X-ray radiations fields 
but low gas densities when compared to PDRs. This in- 
dicates a more efficient gas heating in XDRs than PDRs. 
The high kinetic temperatures resulting from our LVG 
analysis for the AGN sources fit nicely into this picture 
and are even more supported by recent results found in 
nearby AGN through the CO(J=3 2) emission. Mat- 
sushita et al. (2005) detect stronger central emission 
through CO(J=3-2) than through CO(J=1-0) in their 
SMA sample of several nearby Seyfert galaxies. This may 
indicate that temperature effects play a non-negligible 
role for the high HCN-to-CO(J=l-0) intensity ratios in 
AGN. This would also explain the decrease of the HCN- 
to-CO intensity ratios with increasing rotational number 
J in our sample for the AGN sources. However, we also 
find significantly different HCN abundances that may 
support an HCN enhancement through the strong X-ray 
radiation field in AGN in addition to a temperature ef- 
fect. 

6.3. Implication for the SFR - dense gas relation in 
galaxies 

iGao fc Solomon I (|2004al |bl) find a strong correlation 
between the star formation rate (SFR), being propor- 
tional to the infrared luminosity, and the dense gas mass, 
traced by the HCN luminosity. However, the results on 
the AGN dominated galaxies of our sample may indi- 
cate that this correlation is violated in certain AGN en- 
viron ments (compare also lKohno. Nakanishi. fc Imanishil 
12001 . The HCN-to-CO(J=l-0) ratios for NGC 1068 



NGC 5194 and at the center of NGC 6951 dKrips et al 



l2007h do not fall on the correlation of IGao fc Solomon 
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(|2004alTbh . This may be a consequence of the overabun- 
dance of HCN in the AGN sources of our sample. As 
a potential alternative, the HCN-to-CO ratios at higher 
transitions may be better suited for AGN environments 
than the HCN-to-CO (J= 1-0) ratios as they seem to co- 
incide with the SFR-HCN correlation. 

7. SUMMARY AND CONCLUSIONS 

We observed the centre of 12 nearby active galaxies in 
several transitions of HCN and HCO+ with the IRAM 
30m telescope. The results can be summarised as follows: 

1. ) We find that the HCN intensity ratios vary sig- 

nificantly with activity type, i.e., depending on 
which power source dominates the central emis- 
sion. The HCN (R^ /10 ), HCO+ (R? u C j°/ 10 ) and 

HCO + -to-HCN (R^ C °//o HCN ) intensity ratios seem 
to increase with increasing SB contribution, in 
agreement with predictions from theoretical mod- 
els. The highest intensity ratios are found in 
the evolved PDR-dominated starburst galaxy M82 
while the lowest ratios are found in NGC 1068, a 
pure AGN source with a potential central XDR. 

2. ) We also find a variation in the intensity ratios 

among the starburst sources of our sample which 
may be explained by the evolutionary phase of the 
starburst, i.e., a differing dominance of shocks (pre- 
SB), hot cores (young SB), PDRs and SNe/cosmic 
rays (evolved SB). 

3. ) An LVG analysis of the HCN and HCO+ data 

suggests that the SB dominated sources in our 
sample have high molecular gas densities around 
n(H2)~10 40 -10 6 5 cm -3 , kinetic gas tempera- 
tures of T k w20-120 K and HCN abundances of 
Z(HCN)«(0.001-2)xlO~ 8 , while the AGN domi- 
nated regions seem to show n(H2)<10 45 cm -3 and 
temperatures of Tk>40 K with HCN abundances 
of Z(HCN)«(0.1-10)xlO- 7 . 

4. ) The low HCO + -to-HCN intensity ratios found in 

the AGN sources of our sample seem to make it un- 
likely that non-collisional excitation plays a signifi- 
cant role in AGN for the HCN and HCO + emission. 
This may be further supported by the decreasing 
HCO + -to-HCN intensity ratios with increasing ro- 
tational number J in the AGN sources of our sam- 
ple. 



5. ) Assuming thus purely collisional excitation in 

AGN, we can exclude gas density effects as main 
cause for the higher HCN-to-CO (J= 1-0) ratios 
found in AGN favorising hence an increased HCN 
abundance and/or temperature effects in AGN. 
The latter is supported by decreasing HCN-to-CO 
ratios with increasing rotational number J. 

6. ) To explain the potential differences in gas densities 

between AGN and SB environments of our sam- 
ple, we favor a scenario in which the molecular gas 
may be either significantly less clumpy, i.e., has a 
lower 'clump/GMC number density' in AGN than 
SB environments or, alternatively, is more conti- 
nously smeared over AGN environments, i.e., has 
indeed a lower gas density. 

7. ) An estimate of the dense molecuar mass 

to HCN luminosity ( MH 2 -to-LHCN) 
conversion factor X(HCN) results in 
X(HCN)=10±7° Mq (K km s" 1 pc 2 )" 1 for 
the AGN sources in our sample. This is 
a factor of ~2 lower than tha t found by 
iSolomon. Downes fc Radford I (| 19921) : iGreve et all 
(120071) for ULIRGs and starburst sources of 
X(HCN)=20± 3 ° M e ( km s" 1 p c 2 )" 1 , but consis- 
tent w ith the value favored by iGao fc Solomon 1 
(2004^) for their sample of nearby active galaxies 
and ULIRGs. 

8. ) The overabundance of HCN found in the AGN 

sources of our sample indicates that the correlation 
between SFR and HCN luminosity may be violated 
in the vicinity of an AGN. Given the decreasing na- 
ture of the HCN-to-CO ratio with increasing tran- 
sition, we suggest to use the J=3-2 transition of 
HCN and CO as an alternative. 



The results in this paper are based on observations car- 
ried out with the IRAM 30m telescope in Spain. IRAM 
is supported by INSU/CNRS (France), MPG (Germany) 
and IGN (Spain). We are grateful to the dedicated staff 
at the IRAM 30m telescope who have been very sup- 
portive and helpful during the observations. We thank 
the anonymous referee for very helpful suggestions on 
improving the paper. 
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TABLE 1 
Basic properties of the sample. 



Source 


z 


RA 


Dec 


Type- 




NGC1068 
NGC5194 


0.0038 
0.0015 


02:42:40.7 
13:29:52.7 


-00:00:47.9 
+47:11:42.6 


AGN 
AGN 




NGC4826 
NGC3627 
NGC4569 
NGC6951 


0.0014 
0.0024 
0.0008 
0.0047 


12:56:43.8 
11:20:15.0 
12:36:49.8 
20:37:14.5 


+21:40:58.9 
+12:59:29.5 
+13:09:46.3 
+66:06:19.7 


AGN+SB? 
AGN+SB 
AGN+SB 
AGN+SB 




NGC6946 
NGC2146 
M82 


0.0002 
0.0031 
0.0007 


20:34:52.3 
06:18:37.7 
09:55:52.2 


+60:09:14.2 
+78:21:25.3 
+69:40:46.9 


SB 
SB 
SB 




NGC6240 
Mrk231 
Arp220 


0.0245 
0.0422 
0.0181 


16:52:58.9 
12:56:14.2 
15:34:57.1 


+02:24:03.4 
+56:52:25.2 
+23:30:11.5 


ULIRG(AGN+ 
ULIRG(AGN+ 
ULIRG(AGN+ 


SB) 
■SB) 
■SB) 



a type of activity found in the central ±15". AGN denotes 
thereby Seyfert or LINER galaxies. Please note that NGC 1068 
for instance also contains a SB ring but on much larger scales 
(i.e., ~ ±20"). NGC 6951 also houses a strong SB which is still 
within the 30m beam as opposed to NGC 1068. Most clas- 
sifications a re from the NED; com posite nature (AGN±SB): 
NGC 3627 (|Dahlem et all |1994(@D : NGC 4569 (Gabel et al. 
2002); NGC 4826 (Garda-Burillo et al. 2003); NGC 6951 (see 
text); NGC 6240 (Tecza et al. 2000; Gallimore & Beswick 2004); 
Mrk231 (Smith et al. 1999). 



TABLE 2 

Observational parameters as 
taken from the iram 30m 
homepage. 



Observed 


F cf f 


B cf f 


^bcam 


Frequency 






(") 


89 GHz 


0.95 


0.78 


29.5 


177 GHz 


0.93 


0.65 


14.0 


267 GHz 


0.88 


0.46 


9.5 
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TABLE 3 

Line parameteres for the HCN and HCO+ data. 



Source 


tHCN a 

K kms" 1 


tHCN a 

K kms -1 


tHCN a 
x 32 

K kms" 1 


Av b 
km s — 1 


el 0c 

["] 


e 2 s lc 
[") 


["] 


fio d 


f21 d 


f 3 2 d 


NGC1068 


24.5±0.9 

4. / ±U.Z 


20.0±0.4 


19.0±0.6 

<Z.U 


220+10 
i on_l_ 1 n 


4.5 

1 tr, 

lo 


3.0 

i n 
1U 


3.0 

1U 


0.09 
U.Do 


0.16 
U.O / 


0.29 

U.oo 


NGC4826 
NGC3627 
NGC4569 
NGC6951 


6.0±0.4 
2.7±0.2 
2.8±0.1 
3.1±0.1 


4.7±0.4 

<2.0 
2.6+.0.4 
5.0±0.7 


3.4+0.4 

<3.0 
2.3+0.5 
3.6+1.0 


300±20 
290+30 
210+30 
300+20 


20 
8 
10 
17 


20 
8 
10 
17 


20 
8 
10 
17 


0.67 
0.23-f 
0.32-^ 

0.59 


0.91 
0.56 / 
0.67-^ 

0.83 


0.91 
0.71 f 
0.83^ 

0.91 


NGC2146 
M82 c 


q Q_|_n 1 

5.0±0.1 
29±0.2 


1 1 Q-|_n R 

4.4±0.3 
26±0.6 


Q 9-1-0 (\ 

4.3±0.4 
27±1 


1 ^0-1-^ 
lOUZCO 

290+10 
130±5 


JLU 

20 


i n 

1U 

20 


20 


U.OZ J 

0.67 / 
1 


U.O f 

0.91 / 

1 


U.OO J 

0.91^ 
1 


NGU6240 
Mrkzol 
Arp220 


3.2±0.2 
1.9±0.1 
9.7±0.4 


8.8±0.3 
4.9+1.0 
28.4±0.7 


12.3±0.8 
9.3+0.6 
43. Oil. 


410±20 
zzU±zU 
530+20 


3 
3 
2 


3 
3 
2 


3 
3 
2 


0.07 
0.04 
0.02 


0.24 
0.16 
0.08 


0.42 
0.29 
0.15 




THCO+ a 
L 10 




THCO+ a 
x 32 


ZAV 






SilUK' 


els above 








K kms" 1 




K kms" 1 


K kms" 














NGC1068 
NGC5194 


14.6±0.2 
2.4±0.1 




7.6+0.8 
<1.3 


234+4 
134+7 














NGC4826 
NGC3627 
NGC4569 
NGC6951 


3.5±0.1 
2.7±0.2 
2.3±0.2 
2.2±0.1 




<2.5 
<2.1 
<2.1 
<2.9 


300+30 
230+15 
211+17 
300+20 














NGC6946 
NGC2146 
M82 


8.5±0.1 
6.3±0.2 
40.2±0.2 




8.6±0.6 
5.2±0.6 
34.0±2.0 


150±3 
300+20 
130+10 














NGC6240 
Mrk231 
Arp220 


5.0±0.5 9 
1.6±0.2« 
4.6±0.5 9 




8.0+1. s 
3.8+0.49 
8.8+0.99 

















a Velocity integrated intensities, not yet corrected for filling factors; for a definition see equation 2. Error 
are purely statistical and were determined from the gaussian line fits to the data. b FWHM of HCN and 



HCO+ line. c Most of the source sizes have been determined either from the HCN maps where accesible 
(e.g., NGC 6951: Krips et al. 2007b) or from CO maps (NUGA project: Garcfa-Burillo et al. 2004; 
BIMA-SONG: Heifer et al. 2003). Individual galaxies not found in the BIMA SONG or NUGA survey: 
NGC 2146 (Greve et al. 2006); Mrk231, Arp220 (Downes & Solomon 1998); NGC 6946 (Schinnerer et 
al. 2007); NGC 6240 (Iono et al. 2007). d For a definition see Equation 3. While for the first one the 
major axis (FWHM) of the emission was chosen, the minor axis was taken for the elliptical case. e All 
HCN and HCO + lines were averaged over the same area of 30" so that they all should have the same 
filling factor assuming that they come from the same emission region. However, it is not unprobable that 
the higher-J HCN and HCO + lines originate from a more compact region than the HCN(J=l-0) and 
HCO + (J=1-0) line emission but this would make the ratios listed in Table [4] for M82 even higher/ We 
assumed elliptical source sizes for the filling factor here (see definition in Equation 3). The minor axis 
of the emission was taken for the source size. g HCO + data taken from Gracia-Carpio et al. 2006. We 
assumed a 10% error for the integrated intensities here. 
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TABLE 4 
Intensity and luminosity ratios. 



Source 


dHCN a 

■"■21/10 


pHCN a 
-"-32/10 


R HCO+a 
-"32/10 


R HCO+/HCN a 

-"10/10 


„HCO+/HCN a 

"■32/32 


NGC1068 


0.44±0.02 


0.21±0.01 


0.14±0.02 


0.60±0.01 


0.38±0.07 


NGC5194 


0.3±0.1 


<0.23 


<0.2 


0.72±0.06 


<0.7 


NGU4826 


0.38±0.04 


0.23±0.01 


<0.3 


0.59±0.04 


<0.7 


NGC3627 


<0.4 


<0.4 


<0.3 


1.0±0.1 


<0.7 


NGC4569 


0.5±0.1 


0.4±0.2 


<0.4 


0.85±0.07 


<0.9 


NGC6951 


0.7±0.1 


0.4±0.1 


0.44±0.02 


0.70±0.05 


<0.8 


NGC6946 


0.64±0.05 


0.42±0.03 


0.45±0.03 


0.86±0.03 


0.94±0.08 


NGC2146 


0.67±0.04 


0.57±0.05 


0.51±0.06 


1.30±0.08 


1.2±0.2 


M82 (aver) 6 


0.88±0.04 


0.92±0.07 


0.96±0.05 


1.40±0.10 


1.30±0.10 


M82 (ePDR) 6 


1.00±0.06 


1.00±0.10 


1.10±0.10 


1.40±0.10 


1.20±0.30 


M82 (centre) 6 


0.70±0.1 


0.50±0.05 


0.40±0.10 


1.50±0.04 


1.30±0.10 


NGC6240 


0.67±0.05 


0.45±0.04 


0.19±0.07 c 


1.60±0.20 c 


0.70±0.20 c 


Mrk231 


0.62±0.07 


0.58±0.05 


0.30±0.10 c 


0.87±0.09 c 


0.40±0.10 c 


Arp220 


0.69±0.03 


0.50±0.02 


0.21±0.06 c 


0.48±0.05 c 


0.20±0.05 c 


LITERATURE DATA 


NGC253 d 




<0.7 


0.9±0.2 


0.8±0.1 


>1.0 


IC342 e 


0.8±0.05 


0.4±0.04 




0.9±0.1 




NGC4945 / 




0.7±0.1 


0.4±0.1 


l.liO.l 


0.7±0.1 



a corrected for beam filling effects using Equation I (see also Table [3}. b Intensity 



ratios of M82 averaged over the entire disk (=aver), at the position (=(+12", +8")) 
of the eastern PDR (=ePDR) and at the center (=(0",0"))- The ratios for 
the eastern PDR are given because of consistency reason. HCO + values taken 
from Gracia-Carpio et al. 2006. d SB gal axy. Data tak e n from Ma rtin et al. in 
(prep .?). 6 SB galaxy. Data t aken from iNguyen et ahl (pL994) and ISchulz et al.l 
teoOll V SB galaxy, taken from lWang et al.l (|200^1 
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TABLE 5 
Intensity and luminosity ratios. 



bource 


dCO a 
""-21/10 


^HCNIO b 

-kcoio 


^HCN21 b 
-^0021 


^HCN32 b 
£CC>32 


NGC1068 


0.4 


1.9 


1.1 


0.5 


NGC5194 


0.6 


0.5 


0.2 


<0.3 


NGC4826 


0.5 


0.4 


0.3 




NGC3627 


0.6 


0.1 


<0.1 


<0.1 


NGC4569 


0.6 


0.2 


0.1 




NGC6951 


0.7 


0.2 


0.2 




NGC6946 


0.4 


0.2 


0.3 


0.3 


NGC2146 


0.3 


0.01 


0.02 


0.03 


M82 


0.9 


0.2 


0.2 


0.1 


NGC6240 


0.3 


0.1 


0.1 




Mrk231 


0.1 


0.1 


0.3 




Arp220 


0.2 


0.1 


0.3 


0.3 



a corrected for beam filling effects using Equation 1 
(see also Table [3}. b CO luminosities were taken from 
IRAM 30m observations of NGC 3627, NGC 4569, 
NGC 4826 and NGC 6951 in context of the NUGA 
project. The CO(1-0) and CO(2-l) fluxes for 
NGC 1068 were derived from interferometric maps 
(taken from Schinnerer et al. 2000 and Krips et al., 
2007a) while the CO (3-2) emission has been taken 
with APEX during its commissioning which is con- 
sistent with published values. The rest has been 
taken from the literature: Arp220 and Mrk231 from 
Radford et al. 1991 (CO(1-0) and CO(2-l)) and 
Narayanan et al. 2005 (CO (3-2)); NGC 6240 from 
Solomon et al. 1997 (CO (1-0)) and Tacconi et al. 
1999 (CO(2-l)); NGC 2146 and NGC 6946 from 
Gao & Solomon 2004a (CO(1-0)), Braine et al. 1993 
(CO(2-l)), Casoli et al. 1991 (CO(2-l)) and Dumke 
et al. 2001 (CO(3-2)); NGC 5149 from Nakai et al. 
1994 (CO(1-0) and CO(2-l)) and Dumke et al. 2001 
(CO(3-2)); M82 from Walter et al. 2002 (CO(1-0)), 
Thuma et al. 2000, (CO(2-l)) and Dumke et al. 1994 
(CO(3-2)). CO(3-2) data on NGC 3627 are from 
Petitpas et al. (2005). 



TABLE 6 

Results of the LVG analysis in dependency of the dominant activity type. 



Example 




n 


(H 


2) b 


N(HCN)/(dv) c 


Z(HCN)/(dv/dr) d 


[HCN]/[HCO+] c 


source 


K 


cm 


-3 


cm 2 


km s 


km 


s pc 




NGC1068 


20 


10 4 


.o_ 


TO 4 ' 5 


10 15 


-10 16 


10" 


-8. 


-10- 7 


50 




60-240 


10 2 


o_ 


1Q4.S 


10 13 


-10 le 


10" 


-8. 


-10~ 6 


10 


NGC6951 


20 


10 4 


,i> 


-io°- u 


10 lb 


-10" 


10" 


-8. 


-10-'' 


50 




20-120 


10 4 


o_ 


TO 6 - 


10 13 


-10 le 


10" 


-9. 


-lo- 7 


10 


M82(aver) 


60-100 


10 b 


o_ 


10 Bb 


10 ib 


-10 16 


io- 


10 


-10-" 


0.01-1 


NGC6240 


20-120 


10 3 


. >:> 


-10 4b 


10 ib 


-W Lti 


10" 


-7. 


-io~ b 


10 



a Kinetic gas temperature. 1 " H2 density. c HCN column density. d HCN abundance per 



velocity gradient. Abundance ratio between HCN and HCO 



14 



Krips et al. 



- HCNI1-0) 
'-- HCN(2-1) 
HCN(3-2) 




Av [km/5] 



,-0.01 



1- HCN(I-O) 


NGC3627- 


•--HCN(2-1) i~! 




; HCN(3-2) i | 




■ 1 





"- hcn(i-o) 


NGC5194; 


I-- HCN(2-1) 




- HCN(3-2) 


if 1 1 : 




. . 


-500 


500 




Av [km/s] 




Av [km/s] 



Av [km/s] 



_ - HCN(1 - 


-0) '| 


NGC6946 _ 




-- HCN(2- 








HCN(3- 






0.02 




■7 


Li.' 


s? 








"0.01 








1 










M 


1 f \ ft 









'- HGN(1-0) 


NGC6951 \ 


--- HCN(2-1) 




CK(3 2) 




1 . 






:.. |k | 





Av [km/s] 



-500 

Av [km/s] 



Av [km/s] 



Av [km/s] 



-500 


500 


Av [km/s] 


--HCNO-0) 


MB2j 


--- HCN(2-1) 




; ■ HCN(3-2) 









-HCN(1-0) 

- -HCN(2-1) 

- HCN(3-2) 


Mr<231 


0.1 


'-HCNO-0) 
HCN(2-1) 
~-HCN(3-2) I 




' Arp'220; 


. 1 


x ;nrMl^.iir\ 


"0.05 

1 




g 

, , , * i! , , , 




, ii ' , , , 


-500 


500 




-500 





500 



Fig. 1. — Line spectra of HCN(J=l-0) (solid black), HCN(J=2— 1) (dashed red; please note that color figures are only available in the 
online version of the paper) and HCN(J=3-2) (dotted blue) obtained at the centre of each of the 12 galaxies with the IRAM 30m telescope. 
The velocity scale is relative to the LSR velocity of the respective galaxy. Temperature scale is in main beam temperature (Kelvin) and 
has not been corrected for beam filling factors. 
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Fig. 2. — Line spectra of HCO+(J=l— 0) (solid black) and HCO + (J=3-2) (dotted blue) obtained at the centre in 9 galaxies of our sample 
with the IRAM 30m telescope. The velocity scale corresponds to the LSR velocity of the respective galaxy. Flux scale is in main beam 
temperature (Kelvin) and has not been corrected for beam filling factors. The HCN and HCO+ intensities of NGC 6240, Mrk231 and 
Arp220 were taken from Gracia-Carpio et al. (2006) and their spectra are thus not shown here. 
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Fig. 3.— HCN(J=2-l)-to-HCN(J=l-0) and HCN(J=3-2)-to-HCN(J=l-0) intensity ratios of all 12 sources (this paper, Table E} and 
IC 342 (filled grey box; data taken from the literature, Table |S). NGC 1068 (filled blue circle), NGC 6951 (filled green triangle), Arp220 
(open green star), NGC 6240 (filled yellow star) and M82 (filled red box) are highlighted. The dotted grey lines should guide the readers 
eyes and indicate the putative different locations of the SB and AGN dominated sources in this diagram. 
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Fig. 4.— HCO+(J=3-2)-to-HCO+(J=l-0) and HCN(J=3-2)-to-HCN(J=l-0) intensity ratios of all 12 sources (this paper, TableHJl and 
NGC 253 and NGC 4945 (filled grey box; taken from the literature, Table g). NGC 1068 (filled blue circle), NGC 6951 (filled green triangle), 
Arp220 (open green triangle), NGC 6240 (filled yellow star) and M82 (filled red box) are highlighted. The dotted grey lines should guide 
the readers eyes and indicate the putative different locations of the SB and AGN dominated sources in this diagram. 
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Fig. 5.— HCO+-to-HCN(J=l-0) and HCO+-to-HCN(J=3-2) intensity ratios of all 12 sources (this paper, Table!} and NGC 253 and 
NGC 4569 (filled grey box; taken from the literature, Table gj. NGC 1068 (filled blue circle), NGC 6951 (filled green triangle), Arp220 
(open green triangle), NGC 6240 (filled yellow star) and M82 (filled red box) are highlighted. 
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Fig. 6. — x 2 -test results of the LVG analysis (see text for explanation) carried out with MIRIAD including the HCN and HCO + data 
for different kinetic temperatures and abundance ratios of [HCN]/[HCO + ]=0. 01, 0.1,1 (see also Fig. [7]for the higher abundance ratios of 
10 & 50). The dotted red, solid blue, solid green and solid yellow contours correspond to M82, NGC 1068, NGC 6951 and NGC 6240 
respectively. The area enclosed by the contours represent the lowest (reduced) x 2 -test with x 2 ^$l> f° r these low [HCN]/[HCO + ] abundance 
ratios (i.e., <1) we only find a significant result for the SB dominated sources in our sample. The grey lines indicate Z(HCN)/(dv/dr) with 
10 — 5 , 10 -6 , 10 -7 , 10~ 8 and 10 -9 km -1 s pc (from top to bottom, decreasing in line thickness). 
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Fig. 7.— Same as Fig.[6]but for abundance ratios of [HCN]/[HCO+]=10 & 50. 



